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Abstract: This article is a detailed version of our short letter ‘On equilibrium partition 
function for non-relativistic fluid’ [1] extended to include an anomalous ?7(1) symmetry. We 
construct a relativistic system, which we call null fluid and show that it is in one-to-one 
correspondence with a Galilean fluid living in one lower dimension. The correspondence is 
based on light cone reduction, which is known to reduce the Poincare symmetry of a theory to 
Galilean in one lower dimension. We show that the proposed null fluid and the corresponding 
Galilean fluid have exactly same symmetries, thermodynamics, constitutive relations, and 
equilibrium partition to all orders in derivative expansion. We also devise a mechanism to 
introduce ?7(1) anomaly in even dimensional Galilean theories using light cone reduction, and 
study its effect on the constitutive relations of a Galilean Fluid. 
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1 Introduction and Summary 


Non-relativistic systems enjoy an active interest in the physics community primarily for two 
main reasons. First, they are expected to be realised in the low energy physics experiments. 
Second and more fundamentally grounded reason is that a non-relativistic system can be 
thought of as an effective low energy description of an underlying relativistic theory. Hence, 
it is natural to expect that the constitutive relations of a non-relativistic fluid, obtained as an 
effective description of a relativistic theory, may contain new terms which are not considered 
in the coarse grained description of hydrodynamics [3]. For example, if the system breaks 
parity symmetry at the microscopic level, that will enforce us to add parity-odd terms in 
the constitutive relations. The goal of this paper is to revisit the paradigm of non-relativistic 
charged hydrodynamics. We devise a consistent mechanism to derive the parity-even and odd 
terms in the constitutive relations of a non-relativistic fluid up to leading order in derivatives, 
starting from a relativistic theory. 

Galilean fluids^ have been an interesting and active topic of research in recent years [2, 4-11]. 
[2] worked out a consistent way to write Galilean fluid constitutive relations in Newton- 
Gartan covariant formalism, and used the second law of thermodynamics to constrain the 
hydrodynamic transport in two spatial dimensions. It is known that Newton-Gartan geometry 
with Galilean isometry follows from light cone reduction of a relativistic geometry in one 
higher dimension [12, 13]. The idea behind this is that the Poincare algebra in (d + 2)-dim 
has a {d + l)-dim Galilean subalgebra embedded into it. As suggested in [6], this approach 
can be used to construct Galilean covariant tensors in Newton-Gartan formalism, which is 
otherwise a non-trivial task. Similar ideas were also used in [2] where authors constructed 
an extended representation of Galilean group by embedding it in one higher dimension, and 
used it to present the Galilean fluid dynamics in a manifestly covariant manner. 

In this work, we take this approach a step ahead and ask if we can construct a relativistic 
fluid, whose symmetry algebra when restricted to the Galilean subalgebra, is equivalent to 
a Galilean fluid in one lower dimension. This idea has been explored in the past, starting 
with [14] which showed that dynamics of a relativistic fluid reduces to that of a Galilean 
fluid under light cone reduction. However in [15], we observed that this naive approach runs 
into some troubles - the thermodynamics that the reduced Galilean fluid follows is restricted 
(mass chemical potential is not an independent variable; look footnote 7). We also found 
that the parity-odd sector only survives if the fluid is incompressible and is kept in a constant 
magnetic field. It strongly hints that to get the most generic Galilean fluid via light cone 
reduction, we need to start with a modihed relativistic system. 

non-relativistic system is defined by c —^ oo limit of a relativistic system, while a Galilean system is one 
whose isometry group is Galilean. The two are only approximately the same. In this paper however, we only 
talk about Galilean theories, as they are much easier to handle. 
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More precisely, we start with a flat background (metric and gauge field), 

d 

dsjiat = -2dx~dt + ^(dx*)^ Aflat = 0, (1.1) 

i=l 

which has {d + 2)-dim Poincare invariance, (d + l)-dim Galilean algebra sits inside Poincare 
~ all generators which commute with P_ = d- (c.f. [4]). Hence a theory on this background 
which respects x~ independent isometries x^ x^ + ^^{t,x), x^ = {x~,t,x^} enjoys 
Galilean invariance. Gompactifying the x~ direction, we can recover a (d+l)-dim flat Galilean 
background on which non-relativistic theories can be defined. This is known as light cone 
reduction (LCR). We turn on x~ independent fluctuations around the flat background, 

ds^ = —2e“^(dt + ajdx*)(dx“ — Btdt — Hjdx*) + gtjdx^dx^, A = Atdt + Aidx^. (1.2) 


Galilean theories can then be described by a partition function Z[Bt,BiA^O'i^9ijjAt,Ai]. 
Treating these fluctuations as sources, we can define the following observables, evaluated in 
absence of sources (i.e. on flat background). 


P = 


5W 

Wt 


= 

do 


flat 


sw 


€ = 


flat 


5W 

5 $ 


jl = 


flat 


sw 

Scii 


flat ^9ij 


flat 


sw 


q = 


SAt 


flat 


Jq 


sw 

SA. 


flat 


(1.3) 


Here W = InZ, and p, jp, e, j*, A, q, are mass density, mass current, energy density, 
energy current, stress tensor, charge density and charge current respectively of the Galilean 
theory. Invariance of partition function under x~ independent diffeomorphisms will imply the 
following conservation equations. 


dtP + dijp = 0, dt€ + dijl = 0, dtfp + djA = 0, dtq + dijl = 0. (1.4) 

These are exactly what we expect for a Galilean system, if we identify t with Galilean time, 
as suggested by the notation. 

Above procedure can be made manifestly covariant as proposed in [12] and later developed by 
[13, 16, 17] and many others. Consider a curved background, a metric Gmn and a gauge field 
Am with a null Killing vector normalized as V^Vm = 0 and Am = 0; we call this 
background null background‘d. Theories on null background, which we call null theories, are 
demanded to be invariant under preserving diffeomorphisms and gauge transformations. 
The background given in eqn. (1.2) with respective Galilean symmetry then just follows by 
a choice of basis x^ = such that V = d- {t is not necessarily null). It suggests 

^Backgrounds admitting a covariantly constant null Killing vector are termed as Bargmann structures 
in [12], Null backgrounds are special cases of Bargmann structures on which gauge fields and connections’ 
component along the null Killing vector are fixed. 
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that null theories are entirely equivalent to Galilean theories, and are related by merely this 
choice of basis. More formally, null theories exhibit Galilean invariance upon null reduction, 
i.e. getting rid of V direction through compactification. 

We can now study a null fluid on this null background, with the hope to get the most generic 
Galilean fluid after reduction. Unlike ‘usual’ relativistic fluids, in this case isometry V is also 
a background field and hence must be considered while writing the respective constitutive 
relations. This simple consideration happens to resolve all the issues we enlisted before. 
In fact it does much more that that; even before LGR, {d + 2)-dim null fluid is essentially 
equivalent to a (d + l)-dim Galilean fluid, as they have same symmetries. As we shall show, 
their constitutive relations, conservation equations, thermodynamics etc. match exactly to 
all orders in the derivative expansion. 

Another motivation to study null backgrounds is Galilean anomalies^. The thumb rule for 
anomalies tells us that they can only exist in even dimensions. But since light cone reduction 
reduces dimension of the theory by one, even if we start with an even-dimensional anomalous 
relativistic theory, the reduced Galilean theory is odd-dimensional and hence all anomalous 
terms should vanish‘d. However this argument about dimensionality can be bypassed by 
working in null backgrounds. Since there is an extra vector field V in the theory, the tensor 
structure allows for anomalies only in odd dimensions. In fact, one can reconstruct anomalies 
in an even dimensional Galilean theory by starting with an odd dimensional anomalous null 
theory. We would like to use the null fluid construction to see how these anomalies affect the 
Galilean hydrodynamic transport. 

It is known that constitutive relations of a relativistic fluid at local thermodynamic equilibrium 
can be obtained from an equilibrium partition function up to some undetermined ‘transport 
coefficients’ [20, 21]. These coefficients can be determined either from experiments or through 
a microscopic calculation. If we think of Galilean fluid as a limit of an underlying relativistic 
theory, we would expect that its constitutive relations will also follow from such an equilibrium 
partition function, which has been discussed in [7]. We expect that a similar partition function 
can also be achieved via light cone reduction by setting the theory on background eqn. (1.2) 
to be independent of t direction. In this configuration symmetries of the theory break down 
to diff X U{1)^ (spatial diffeomorphisms, Kaluza-Klein transformations, mass transformation 
and gauge transformation), and one can easily write down the equilibrium partition function 

®We will only be talking about global t’Hooft anomalies appearing in Galilean theories as described by 
[18]. Our working definition of anomaly shall be that the respective conservation laws are violated by certain 
terms purely dependent on the background sources. We do not dwell in the microscopic interpretation of these 
anomalies. 

"^This is in contrast with the results of [18], where author found that the relativistic anomalies survive the 
light cone reduction and show up as gauge/gravitational and Milne anomalies in the Galilean theory. We 
observe that this is because of the presence of extra scalar sources in the Galilean theory that are reminiscent 
of reduction and must be switched off in a physically realizable theory. We present a detailed analysis on these 
issues in a companion paper [19]. 
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invariant under these symmetries as a gauge invariant scalar made out of the background 
fields. In null background picture, same story follows by introducing another isometry , 
and choosing a basis such that K = dt- 

Hence the refined goal of this paper is to set up a consistent theory of hydrodynamics on 
null backgrounds. We want to find the most generic constitutive relations for a fluid on 
null backgrounds constrained by the second law of thermodynamics and requirement of an 
equilibrium partition function. Later employing light cone reduction, we interpret these null 
fluid constitutive relations as constitutive relations of the most generic Galilean fluid. 

Organization of this paper is as follows. In § 2 we review the construction of torsionless null 
backgrounds, and construct an equilibrium partition function for null theories. Then in § 3 
we study hydrodynamics on these null backgrounds, and put constraints on its dynamics by 
equilibrium partition function and second law of themodynamics. We devote § 4 to review 
the procedure to obtain Galilean theories from null theories by light cone reduction, and use 
it to study Galilean hydrodynamics in § 5. Finally in § 6 we extend this entire construction to 
anomalous fluids. In appendix (A) we extend the entropy current calculation in presence of 
minimal compatible torsion, which is required to get agreement between equilibrium partition 
function and entropy current constraints. In appendix (B) we express all these results in 
conventional non-covariant basis. In appendix (C) we provide a comparison of our results 
with those of [2]. At the end, in appendix (D) we mention notations and conventions of 
differential forms used throughout this paper. 

2 Construction of Null Backgrounds 

We start our discussion by formally setting up null backgrounds, which will prove to be a 
natural ‘embedding’ of Galilean (Newton-Cartan) backgrounds into a spacetime of one higher 
dimension. These kind of backgrounds were first considered in [12] and further explored by 
[13, 16, 17] where authors recovered Newton-Gartan gravity by light cone reduction of general 
relativity. We will refine the approach by constraining the background field content so that it 
exactly matches that of a non-relativistic theory, hence letting us study physically realizable 
Galilean fluids later. 

Let us consider a manifold M-(^d+ 2 ) equipped with a metric ds^ = Cmn^x^ dx^ and a U{1) 
gauge field A. = AMdx^ together referred as background fields/sources. AI(rf_|_ 2 ) is also 
provided with the Levi-Givita connection, 

= 2 ^^^ {QmGns + QsGnm - QnGms) , ( 2 . 1 ) 

and a covariant derivative Vm associated with T^^g and Am- We demand that physical theo¬ 
ries on A4(rf_|_2) are left invariant by diffeomorphisms and gauge transformations parametrized 
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by infinitesimal parameters = {^ = which we call symmetry data. Action of 

'ijj^ (denoted by 6^) on various background fields is given as, 

S^Gmn = £^Gmn = 2V(j\^^jv)5 = 5 m (A(^) + ^^An) + (2.2) 

where denotes Lie derivative along ^ and J^mn is the field strength of Am- One can check 
that symmetry data form an algebra with commutator defined by, 

. (2.3) 


Correspondingly their action on a general field (suppressing all the indices) (p also forms an 
algebra with commutator given by [(^^i, Physical theories on A4(^d+2) can 

be described by a generating functional^ W[Gmn,Am] which is seen as a functional of the 
background sources. Under infinitesimal variation of these sources linear response of W is 
captured by, 


6W 




]^T^^6Gmn + J^^5Am 


(2.4) 


rpMN ^ jM energy-momentum tensor/current and charge current respectively. De¬ 

manding partition function to be invariant under the action of "0^ given in eqn. (2.2), we can 
obtain a set of Ward identities these currents must follow. 


VmT^^ = VmJ^ = 0. 


(2.5) 


These are the energy-momentum and charge conservation laws of a relativistic theory. It is 
not mandatory for a physical theory to admit a Lagrangian description, in which case the 
theory itself can be characterized in terms of conserved currents with dynamics 

provided by equations of motion (2.5). 


2.1 Compatible Null Isometry 

So far whatever we have said applies to any relativistic theory. We now specialize to our 
case of interest - ‘null backgrounds’ by introducing a null Killing vector. More formally, a 
symmetry data ifv = {V = U^5m, A(y)| will be said to generate a compatible null isometry 
on A4(rf+2) h h follows, 

1. Action of 'ifv is an isometry, 5vGmn = SvAm = 0, 

2. V is null, V^Vm = 0, 

3. V is preserved under covariant transport, VmV^ = 0, and, 

4. Component of gauge field A along V is fixed to: V^^Am = ~^{V)- 
® Actually W = In Z where Z is the QFT generating functional. 
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We will call backgrounds admitting a compatible null isometry to be null backgrounds. Since 
we are working with torsionless manifolds, one can check that above conditions imply that 
T^mn = QaiVn — djsfVM = 0. This is a dynamic constraint and can be violated by quantum 
fluctuations off-shell - a fact that will become important when we write equilibrium partition 
function for fluids on null backgrounds in § 2.2.1. 

Null backgrounds possess some nice features, first one being: m‘P = for any contra- 

co-variant tensor ip (all indices suppressed), transforming in appropriate representation of the 
gauge group. Further if ip is entirely made up of by first consistency condition, 

= dvp = 0. These consistency conditions also imply, 

Fnm = V^TZnmrs = TlNMRsy^ = 0 , ( 2 . 6 ) 

where TZmnrs is the Riemann curvature tensor. We term physical theories on null back¬ 
grounds (with compatible null isometry ^jJv) as null theories, and demand them to be invari¬ 
ant under ^jJv preserving symmetry transformations i.e. [^p^,^l;v] = 0. This will break down 
the Poincare symmetry algebra to Galilean, and give null theories a Galilean interpretation. 
Algebraic relations (2) and (4) in the definition of compatible null isometry, will imply, 

KV^Vm) = 0 ^ V^V^SGmn = -2Vm6V^, 

5{V^ Am + \v)) = 0 => V^6 Am = -SA(v)-^m6V^. (2.7) 

It immediately follows that under a variation of background sources restricted by dtpv = Oj 
linear response of partition function eqn. (2.4) is still completely characterized by , 

with an added ambiguity in currents, 

tMN ^ + OiV^V^, ^ + d2V^, ( 2 . 8 ) 

where 0’s are some arbitrary scalars. One can check that they leave the conservation equations 
(2.5) invariant, provided dyOi = 5 y 62 = 0. This point onwards whenever we talk about 
variation, we implicitly assume it to follow 6'ijjv = 0. Also for the following analysis, we 
partially fix 'ipy by choosing A(y) = 0 for convenience. 

2.2 Equilibrium 

For our later discussion on equilibrium partition function of hydrodynamics, it will be helpful 
to define a notion of equilibrium on null backgrounds. A system is said to be in equilibrium 
if it admits a time-like isometry generated by 'ipK = {K = Om, A(^x)}^ be. Km < 0 
and 6kGmn = SkAm = 0. Using K we can define a null field, 

1 A'mK" 

{K +B,V ), e =-J^, 
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which is orthonormal to V, i.e. = —I, = 0- We define spatial slice -M-^d) 

as the spacetime transverse to V and with projection operator, 

( 2 . 10 ) 

Using diffeomorphism and gauge invariance of ^6 pick up coordinates = {x“, t, x®} 

such that, 

V’V = {u = 9-, A(y) = 0} , 'i/jK = {K = dt, A(^k) = 0} , (2-11) 

and coordinates x = {x®} span In this basis, background fields can be decomposed 

3 jS, 


~ GjvfAfdx^dx^ = —2e (dt + a^dx®) (dx 
A- = Am^-x^ = Atdt + .Ajdx®. 


Btdt — iJjdx®) + Qijdx^dx^, 

( 2 . 12 ) 


Indices on can be raised and lowered by and its inverse gij. Decomposition of other 
derived fields follow trivially from here, 



» j 

fe^Bt\ 


(-l\ 

Vm = 


0) 

II 

) ^{K)M = 

Bt 


i^-e ^aij 

Y 0 y 


\B,) 



/o 0 0 \ 


( BkB^ -akB^ B^ \ 

P{K)MN — 

0 0 0 

tdMN 

’ AK) - 

—akB^ a^ak —a^ 


o 

O 


y I?® —a® 5®-^ j 


(2.13) 


where we have defined Bi = Bi — aiBt, Ai = Ai — aAt- Under this choice of basis, one can 
check that residual symmetry transformations are parametrized by x dependent symmetry 
data V'c = ^(c)}) which acts on reduced set of background fields as. 


d*, ^^ai — A £^ ai^ ^^Bf — £^Bt^ S^B^ — “b £^Bi^ S^gij — £^gij^ 

5 ^At = £^ At, ^^Ai = £^ Ai + 

where denotes lie derivative with respect to Its trivial to see that ai,Bi,Ai transform 
as U{1) vector gauge fields, ^,Bt,At transform as scalars, and gij transform as rank 2 tensor. 
The response of partition function eqn. (2.4) in equilibrium under infinitesimal variation of 
these sources can be worked out to be. 


/ {dx®}V5 


(Ti_ + T__Bt) ^ [T\ + JAt] 6 ai 


1 




1 


+—T^gij + T__5wo - —TA6Bi - J_6uo - -^JAAi 




1 




(2.14) 











where we have defined: 


■do = '^o = \t3t, Vo = (2.15) 

■u iy 

d = l/(/3ii) where /3 is the radius of the euclidean time t = it and R is the radius of 
compactified x~ . We define a connection on M.(^d) 

1% = (didji + dkgij - digij), (2.16) 

and Vj as its associated covariant derivative. We call the associated Riemann curvature 
tensor Rijki- Note that condition for torsion-less manifolds, Rmn = 0 implies that in equi¬ 
librium, 

fij = dittj - djtti = 0, dido = 0. (2.17) 

Again, these conditions can be violated offshell, which will be important in next subsection 
when we start construction the equilibrium partition functions. 

2.2.1 Constructing Equilibrium Partition Function 

Motivated by applications in hydrodynamics, we want to write the most generic form of 
equilibrium partition function allowed by symmetries arranged in a derivative expansion of the 
background sources. Partition function is generally written as integration of scalar densities. 
While the partition function is itself invariant under symmetries, such statement cannot be 
made for the integrand. In fact, terms can be added to it whose variation is gauge invariant 
only upto some boundary terms. We can hence decompose into®, 

f = ( 2 . 18 ) 

J Vo JMi^d) 

where Phs is a scalar, and is the d-dimensional ^ Chern-Simons’ form. 1^^ is defined such 
that is a ‘polynomial’ made out of field strengths da, dB, dA and curvature 

2-form Wj = Rkdjdx^ A dx^ It is known that written in odd spatial 

dimensions (d = 2n — 1), and upto first non-trivial order in derivatives is given as, 

n n—r 

'Pcf = E (”) E A {dB^^ A (dda)^("-^-^), (2.19) 

r=0 s=0 

where C^r,s) a-re constants and (”) = is the binomial coefficient introduced for later 

convenience. Note that the torsionlessness condition 'H = 0 would imply f = da = 0 on-shell, 

® Usage of subscripts Hg, Hv is motivated from eightfold classification of relativistic transport in [22]. It 
is yet not clear if such a classification is also applicable to null backgrounds, so for us this usage is purely 
notational. 
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but since partition functions are to be written off-shell, we include these terms. From here 
we can find the on-shell variation of ignoring some boundary terms, 


n—1 


= - 


n 


'M 


AiABr 


n—r—1) 


(d) r=0 


A I + C2,(^+i)(5A + C2Xr)^B 


( 2 . 20 ) 

where C'i^(r) = C(^r,n-r-i) and C 2 ^[r) = C'(r,n-r)- From here we can trivially read out the 
contribution of to currents in equilibrium; we will come back to it in § 3.3.1. Coming 
back to VFhs^, it is now just integration of the most generic scalar Phs rnade out of background 
sources arranged in a derivative expansion. At ideal order (no derivatives) Png^ideai = Po is 
defined as a gauge invariant function of 'do-, 'cuo, ^o-, 


W, 


eqb 


Hg,ideal 


= j {dx*} 


( 2 . 21 ) 


We take this opportunity to define near equilibrium thermodynamics on null backgrounds. 
It is known that euclidean partition function Z = e'^ can be identified with grand canonical 
partition function of statistical mechanics. It follows that P can be identified as pressure 
de-nsity, id as temperature, and w as chemical potentials scaled with temperature, which 
boil down to their o values at equilibrium. Differential of P{'d, w, v) can be expanded as. 


dP = S'di? -|- Rd{wd) + Q<i{v'd), 


( 2 . 22 ) 


where we identify S as entropy density, R as ^mass density^ and Q as charge density. We can 
also define an energy density E by invoking Gibbs-Duhem relation, 


E = Sd + dRw + dQu - P. (2.23) 

Taking a derivative of this relation and using eqn. (2.22) we can find the first law of thermo¬ 
dynamics, 

dE = MS PdwdR + dudQ. (2.24) 

Existence of a mass density makes this thermodynamic system already look Galilean and we 
take it as first hint that (at least) near equilibrium theories on null backgrounds are secretly 
Galilean'. Coming back to equilibrium, we can vary partition function eqn. (2.21) and use 

^ It is interesting to see that a null theory satisfies different thermodynamics than a relativistic theory. It was 
noted in [15] that if we start with a relativistic fluid (following relativistic thermodynamics), thermodynamics of 
Galilean fluid after reduction gets restricted. In our setting this restriction manifests itself as E + P + R'd'CJ = 0. 
After a non-trivial redefinition of thermodynamic functions, 

Erel — ‘2E P, Prel — P, Srel — 5*, Qrel — Q, 'ttrel — ut}, /Trei — where CL — . =, 

a a V—2i?ci7 

this restricted thermodynamics is equivalent to relativistic thermodynamics, 

dEj-^l — aj-^ldS-rel -f pveldQrel, Ej-^l — Srel'^rel H” Qrell-^rel Prel- (2.25) 
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eqn. (2.14) to read out components of currents in equilibrium at ideal order, 


)o,ideal = , (T ) o,ideal = Ro, (Tf- + T Bf) o,ideal = Eg, -{J-)o,ideal = Qo, 

(2.26) 

and rest all spatial currents zero. Using , ^{K)' recompile these into covariant 

language, 


rpMN 
^ o,ideal 




+ 


+ PoPf 


MN 
{K) > 


jM 

^ o,ideal 


QoV( 


M 

{Ky 


(2.27) 


These look like some sort of ideal fluid constitutive relations, but are quite different from a 
relativistic fluid. We will make notion of this fluid on null backgrounds ~ null fluids more 
precise in next section. 


3 Hydrodynamics on Nnll Backgronnds 

Having already developed some intuition in last section, we proceed to formally construct 
hydrodynamics on null backgrounds in its full generality. Any quantum field theory in near 
equilibrium regime can be described by hydrodynamics. Systems having a hydrodynamic 
description, fluids are in local thermodynamic equilibrium, i.e. variations away from equi¬ 
librium are on scales much much larger than the characteristic scale of the system. This 
essentially means that variables defining the fluid are much much larger than their space-time 
derivatives. One can therefore express observables (currents) of the theory as a derivative 
expansion of symmetry covariant data made out of fluid variables. 

Note that conservation laws (2.5) are {d + 3) independent equations, so any system with 
{d + 3) variables would be exactly solvable. We choose to describe our system by a fluid 
with null velocity u normalized as u^um = 0,u^Vm = —1, which will give us (d) degrees of 
freedom, and three thermodynamic variables: temperature d, mass chemical potential •dw, 
and charge chemical potential 'dv. We are interested in configurations which respect the 
isometry generated by V’y- Hence, all the fluid variables as well as constitutive relations are 
annihilated by action of tpv 

Hydrodynamics (due to dissipation) is not described by a partition function; rather it is 
characterized by the most generic form of currents in terms of background fields 

and fluid variables ,'d,w,u known as ‘‘constitutive relations’ of the fluid. 
Dynamics of these currents is given by Ward identities eqn. (2.5) imposed as equations of mo¬ 
tion. The constitutive relations are further constrained by certain physicality arguments, like 

Interestingly this map between relativistic and restricted null thermodynamics is exactly the same as the map 
between relativistic and restricted Galilean thermodynamics found in [15] by null reduction with a = “u”*"” in 
their language. It motivates us to propose that thermodynamic systems on null backgrounds are equivalent to 
thermodynamic systems on Galilean backgrounds. 


- 11 - 



second law of thermodynamics or existence of an equilibrium configuration. Using 
and we can decompose constitutive relations as, 

pMN = TZu^u^ + 

J^ = Qu^ + I^, (3.1) 

where we have used redefinitions eqn. (2.8) to get rid of some terms. P,£,P,Q are some 
arbitrary functions of The tensors E^, T^'^, contains derivative corrections 

and are transverse to and , and is traceless. Comparing constitutive relations 

eqn. (3.1) to eqn. (2.27) we can refer that, at ideal order P,£,P, Q boil down to respective 
thermodynamic variables R,E,P,Q. In the presence of dissipation however, these functions 
can deviate from their thermodynamic values. 

3.1 Hydrodynamic Frames 

Note that fluid variables , i}, w, v are some arbitrary dynamical fields introduced to describe 
the near equilibrium quantum system. Like any field theory, these fields can be subjected 
to arbitrary held redehnition, called the hydrodynamic redefinition freedom. Some of this 
freedom is already hxed by the ideal order equilibrium partition function, requiring that 
these helds boil down to V^y'do,Wo,i'o in equilibrium conhguration at ideal order. Away 
from equilibrium however we are free to perturb these variables the way we like as long as 
the mentioned restriction holds, 

^ + 6u^, ^ d + dd, w^w + 6w, + (3.2) 

where the variations are some arbitrary functions of huid variables and background helds, 
subjected to velocity normalization conditions umSu^^ = VmSu^ = 0. Note that near equi¬ 
librium assumption requires these variations should contain at least one derivative. For the 
physics to remain invariant under these transformations, we require that the functional form of 
j’MN^ jM pgj^ain unmodihed. Hence to hrst non-trivial order in derivatives we will get, 

rj^MN ^ (7^ + 5R)u^u^ + 2{£ + 5E)u^^V^^ + {V + 6P)P^^ 

+ 2 + R6u^^^ + 2 (^e(^ + {E + P)6u^^^ 

^ {Q + 6Q)u^+ Q6u^^+S^, (3.3) 

from where we can obtain hydrodynamic frame transformations, 

n^n + 6R, £^£ + 6E, V^V + 6P, Q^Q + 6Q, 

+ ^ + {E + P)6u^. (3.4) 
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Out of these we can construct three hydrodynamic frame invariants, i.e. quantities that do 
not transform under hydrodynamic frame transformations, 


_ r^AIN _j_ pMN 


(v-p)-(e-E)^p-in-R)^p-(a -« Ap 


TM ^ _ 


M QmM 


gM ^ jgiU _ 


R 


M 


E + P 
R 


>M 


(3.5) 


All the physical information about fluid constitutive relations is encoded in these invariants. 
It is sometimes convenient to fix a hydrodynamic frame to be able to talk about the physical 
constitutive relations directly. Most popular choices involve identifying £,TZ, Q with E, R, Q, 
and dumping all the dissipation into V. This fixes the ambiguity in 'd,TU,n. 


For fixing the velocity redefinition, in spirit with the ‘usual’ relativistic fluids we can use, 
‘Eckart Frame’ in which is chosen to be zero, or ‘Landau Frame’ in which is zero. A 
more natural® frame in this case is the ‘Mass Frame’ where is chosen to be zero, which 

aligns velocity along R flow and all dissipation transverse to . We will mainly work in 
mass frame for which constitutive relations are given as. 


pMN ^ ^ 

= Qu^ + (3.6) 


Another helpful frame for our work is to choose all the fluid variables to be equal to their value 
at equilibrium® exactly, not just at ideal order. We call this ‘Equilibrium Frame’. This has 
the advantage that equilibrium partition function naturally gives constitutive relations in this 
frame. To be precise, in equilibrium configuration setting {u^,Td,w,i'} = 'do, ^ 0 , 1 ^ 0 } 

in the constitutive relations eqn. (3.1), and putting them into equilibrium partition function 
variation eqn. (2.14), we can deduce that, 


TZo = 


6W, 


eqb 

(d) 


S 00 o 


= %- 


5W. 


eqb 

id) 


SBi 


Qo = 


sw, 


eqb 

(d) 


6uo 


= %- 


5W. 


eqb 

id) 


hA,; 


= n- 


5W. 


eqb 

(d) 


6^0 


Ej, - Wo^oK - 


bW\ 


eqb 


id) 


SOi 


V + = 2??o- 


5w: 


eqb 


id) 




(3.7) 




Switching back and forth between frames is a non-trivial task, and has to be done order by 
order in derivatives. We shall see in the subsequent sections that different physical aspects of 
our theory of interest are better understood in different frames. We have to switch between 
frames accordingly. 


®Upon reduction this will imply that mass current does not have any dissipation, i.e. we associate the fluid 
velocity with the flow of mass. 

®This frame choice however does not completely fix the hydrodynamic ambiguity. You can still shift fluid 
variables with terms that vanish in equilibrium. However for equilibrium partition function calculations, it is 
good enough. 
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3.2 Entropy Current 


Since hydrodynamics is an effective field theory, we start by writing down all possible ex¬ 
pressions, compatible with symmetry, that can contribute to E^,In addition 
since we are dealing with a thermodynamic system, we must ensure that the second law of 
thermodynamics is satisfied, i.e., there must exist an entropy current whose divergence 
is positive semi-definite, 

> 0. (3.8) 


We can construct the most generic entropy current for the fluid as. 


jM _ jM I 'Y'M 

^s ^s^can s ’ 


jM 

^s,can 


i) V 


(3.9) 


which is just Su^ at ideal order. J^^an is called the canonical entropy current, and is given 
purely in terms of constitutive relations. on the other hand, are arbitrary derivative 

corrections to the entropy current. Note that T^, unlike T^, is not required to be transverse 
to and . Using hrst order equations of motion we can obtain. 


u^OmE = 


jMN 


R {^NRU 


R 


where we have defined. 


{E + P)Q, u^SmR = -R&, u^OmQ = 
- 'ddN'!^) - {E + P)^dN'& + Q {Enru^ - 'ddNv) 


-ge, 

= 0 , 


^MN = duUN — QnUm, 


0 = Vmw^. 


(3.10) 

(3.11) 


Using these, divergence of canonical entropy current can be computed to be, 

^^MJ^lcan) = [EmNU^ " , (3.12) 

which will come in handy later. Note that each term in above expression is product of 
derivatives (called composites). This heavily constraints the form of T^. Its divergence 
must not contain any pure derivative terms (terms which are not composites), otherwise total 
entropy current cannot be ensured positive semi-dehnite. 


In next subsection, we write the most generic constitutive relations of a null fluid upto lead¬ 
ing derivative order^^ in parity-odd and even sectors. We further impose constraints on these 
constitutive relations by imposing second law of thermodynamics and requirement of an equi¬ 
librium partition function independently, and compare the results from both the approaches. 
Readers who are more interested in Galilean fluid results, can skip this computation and di¬ 
rectly proceed to subsection § 3.4 where the final results for null fluid have been summarized. 
These results can be directly used to read off the constitutive relations of a Galilean fluid, 
which has been done in §§ 4 and 5. 

leading order we mean the first derivative corrections appearing in constitutive relations. In parity- 
even sector it happens at one derivative order itself. In parity-odd sector however it depends on the number 
of dimensions - in odd dimensions {d = 2n — 1) the first correction appear at (n — 1) derivative order, and in 
even dimensions (d = 2n) in appears at (n) derivative order. 
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Table 1: Leading Derivative Order Data for Null Fluid 


Data 

Value at Equilibrium 

Parity Even 

Q = Vmu^ 

0 


0, 0, 0 

P^^Bnw, P^^Bnv 

0, vv, vv 

pMN _ tDdNw) 

0 

P^^ {J^NRU^ - ^BnI^) 

0 

aMN = 2P^Rp^^V^RUs) - iP^^Q 

0 

Parity Odd — Odd Dimensions (d = 2n — 1) 

= ★ [F A n A A D^(«-'-i)] ^ 

III 

* 

> 

■i 

> 

> 

1 

1 

Parity Odd — Even Dimensions (d = 2n) 

= * [F A n A A 

Z,(,)|;^Q^*[(dA)^''A(dH)^(-^)] 

if^fBNW, l^fBNR 

0, 

{^NRU^ - '&BnVo) 

0 

iPf {J^NRU^ - ^Bnu) 

0 

,R(M N) 

(r) R 

0 

where, 

* [F A n A A 



3.3 Leading Order Hydrodynamics 

In [23] we discussed in detail the procedure to count various independent data that appear 
in constitutive relations of usual relativistic fluids. This can be easily extended to null fluid. 
However in this work we are only interested in leading order null fluid, so we can write the 
required data by hand without going into the technicalities of [23]. All possible scalars, vectors 
and symmetric traceless tensors made out of background fields and fluid variables has been 
enlisted in table (1); data marked with * can be eliminated by using first order equations of 
motion (3.10). 

Using data in table (1) we can now write the most generic form of leading order constitutive 


- 15 - 
























relations. For parity-even sector we will get, 


- p^^ce, 

£(^) = [XevodN-^ + XeydNV + K^dN'O + CTe {PnRU^ “ I^Qn^)] , 

pf() = P^^ [Xq^dNZU + XqudNV + Kq^AT'!? + CTq {PnRU^ “ '^Onr)] ■ (3.13) 

In parity-odd sector however, in odd number of dimensions {d = 2n — 1) we will get. 


n—1 


n—1 


n 


MN 

(n-1) 


= 0 , 




M 


(n-1) 


E t?y,) = (”;P,(r|!;r). 


‘r=0 


‘r=0 


(3.14) 


and in even number of dimensions {d = 2n), 

n n—1 

f,MN _ _pMNSr^ (n\^ , iR{^„N) 

^{n) P Z_^ \r)^{r)^r) { r )d(ry(r) ^ R' 

r=0 -r^O 

n—1 

^ X] K'uu(r)dN'^ + Kv{r)dN^ + ^e{r)dN'& + ^e(r) {^NRU 


r=0 

n—1 


R 


fn-l\iMN 

^(n) - Z^[ r )‘(r) 
r=0 


'^Onr) , 

■dd^v) . 

(3.15) 


Similarly we can work out constitutive relations to arbitrary high derivative orders, but in 
this work we will not be interested in those. 


3.3.1 Constraints throngh Equilibrium Partition Function 

The constitutive relations in described above are constrained by the requirement of existence 
of an equilibrium partition function. The statement is that at equilibrium, any null theory 
must be determined by the most generic partition function made out of background fields 
discussed in § 2.2. We have already seen that at ideal order, equilibrium partition function 
gives thermodynamic meaning to various functions. Even at further order in derivatives, equi¬ 
librium partition function turns out to be very useful to (partially) determine the constitutive 
relations. It gives constraints on various transport coefficients, and tells us which of them are 
physical. We along with many people in past have used this approach to find transport of a 
relativistic fluid. Here we attempt to outline a similar procedure for null fluid up to leading 
order in derivatives. 

Leading Order Parity Even Sector: Leading order parity even sector contains one 
derivative corrections to ideal fluid dynamics. Using table (1), we see that at equilibrium 
only terms coupling to A’s survive in frame invariants eqn. (3.13), 

P^o{l) ~ ^o(l) ~ -^oeroV*n7o -|- Aoei/V^t'oj ~ Xoqvj'^^'^o + Aoqi/V^t'o- (3.16) 
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On the other hand there are no one-derivative scalars at equilibrium to construct partition 
function. Hence all the coefficients appearing above must vanish, 


(3.17) 

Since equilibrium partition function is identically zero, none of the fluid variables get order 
one even correction out of equilibrium (in mass frame). 


— ^qvj — — 0 




Leading Order Parity Odd Sector (for d = 2n — 1): In odd dimensions, d = 2n — 1, 
the first parity odd contributions show up at (n — l)-derivative order. At equilibrium all the 
parity-odd terms survive in constitutive relations eqn. (3.14). On the other hand there are no 
gauge invariant scalars to construct equilibrium partition function, and it gets contributions 
only from the Chern-Simons piece (cf. § 2.2.1). Consequently we get constitutive relations in 
equilibrium frame. 


n—l 


®o(n-l) — “ ^ 0 *^ 2 ,(r) “ l^oC2^(r+l)) , 

r=0 
n—l 

r=0 

n—l 


r=0 


(3.18) 


and rest all zero. Here Ci C' 2 ^(r) are constants introduced in § 2.2.1. Performing a hy¬ 
drodynamic frame transformation, we can get the transport coefficients introduced in frame 
invariants eqn. (3.14) as. 


E + P — -dwR 
R 


^e{r) = -I- 

^q{r) • 


C'2,(r) - '&^C‘ 


a.h-ri) 


(3.19) 


We see that both (set of) transport coefficients are completely determined upto some con¬ 
stants. Outside equilibrium, fluid velocity gets a correction (in mass frame) given by: 


Q n 1 


Rn 




o(r)C'2,(r)- 


r=0 


(3.20) 


Corrections to other components of velocity can be determined by this using normalization 
conditions. Other fluid variables do not get any leading order odd correction. 
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Leading Order Parity Odd Sector (for d = 2n): Contrary to the last case studied, in 
even dimensions d = 2n, the first parity odd contributions show up at n-derivative order. 
In even number of dimensions, only terms coupling to A’s and C, survive in frame invariants 
eqn. (3.15), 


n*! . = - 


o(n) 




r=0 


n—1 


d(„, = E ("d)C 


(r) 


r=0 

n—1 




r=0 


h(„) = E ("r'Kr, 


(3.21) 


On the other hand using data in table (1), we can write the equilibrium partition function 
as, 


^ / ^ (-) ^ ("7)5„(,,,)(dA)^^ A (dB)^^ A (Ma)^(-^-^). (3.22) 

dM(d) j.=o s=0 

Varying this partition function on-shell, we can compute the constitutive relations in equilib¬ 
rium frame. 


Po = 0, Tjf' = 0, 

_d 

' di9o 


So = —'&0 ir) TZq = — Y^ ir) 'g^'^o2,{r)lo{r), Qo = — Y^ ir) '^'^o2,{r)^o{r) ^ 


r=0 

n—1 


r=0 


r=0 


K = '^lriY{'^r^)^l\r) {^jSol,{r) “ '^odjS^2,(r) “ ^odjS^2,{T+l)) 

r=Q 

n—1 n—1 

Ri = -tf„«E (”;‘)CwaA.2,(r) Ji = -tfo^E (”;‘)CM«A.2,(r+l), 


(3.23) 


r=0 


r=0 


where 5^ = S(^r n-r- 2 )^ ‘^ 2 ,(r') = ‘5(r,n-r-i)- Transforming these to mass frame, one can 

compute the transport coefficients appearing in eqn. (3.21), 
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We see that 5 (set of) transport coefficients C(r)) -^ero(r)j ^ev{r)^ ■^gro(r)) ^qv(r) determined 
in terms of 2 (set of) functions S 2 ^{^r)- Corrections to fluid variables outside equilibrium 

in mass frame are given as, 


= 

= 

= 

= 


E 

r=0 

n 

E G)Ur) 


' 2 d-do d d-dn d 


+ 


+ 


d-do d 


r=0 

n 


° dEo d-do dRo dwo dQo dvo 

„9 dvuo d dwo d dwo d 
+ — + 


EC)'. 


o(r) 


r=0 


I? 


° dEo d'do dRo dwo dQo dvo 

dvn d dvn d dvo d 


2 '-’>'0 


+ 


+ 


° dEo d'do dRo dwo dQo duo 


^o2,(r )) 


^o2,(r )) 


^o2,(r )) 


A 

Ro 


n—1 




(3.25) 


r=0 


3.3.2 Constraints throngh Entropy Cnrrent 

As we have said, the second law of thermodynamics for null fluid implies the existence of an 
entropy current with non-negative divergence. From our experience of usual relativistic fluids, 
we expect second law requirement to give all the constraints we found through equilibrium 
partition function, and more. However as we shall see, we will not get all the partition 
function constraints through entropy current. This can be accounted to the fact that in 
this computation we will miss constraints coupling to EL = dV, which is set to zero by 
requirement of manifold being torsionless. Since this condition can be violated off-shell, 
equilibrium partition function can however ‘see’ these constraints. In appendix (A) we will 
turn on minimal amount of torsion to allow non-zero EL, and will verify that we get all the 
partition function constraints through entropy current analysis as well. Here we perform 
torsionless computation to leading derivative order. 

Leading Order Parity Even Sector: At leading even order, no terms can be introduced 
in without having pure derivative terms in the divergence, hence, = J^oan) ''^hose 
divergence using eqns. (3.12) and (3.13) is given as, 

{K^dNW + KvdNv) dM'&+P’^^ {K^dNW + KudNv) [Emru^ - '&dMl') 
- i(cr, - l^Kq) [FnRU^ - '&dNv) P^^^du^ 

- ^K^P^^dM'&dN'd + OqP^^ {Pmru^ - '&dMv) {Pnru^ - ’ddNv) + ]^r]a^^aMN + C©^- 

(3.26) 
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Demanding S/mJ^ > 0, from the first line we get all the equilibrium partition function 
constraints eqn. (3.17), and in addition from last two lines, 


C,crg>0, Ke < 0, (Te = 'dKg. 


(3.27) 


Leading Order Parity Odd Sector (for d = 2n — 1): In parity odd sector however, there 
are terms we can write in which have composite divergence. We first consider the odd 
dimensional case for which we will get^^ , 


_ 


n—1 

E 

‘r=0 


n 




i ^Ar ' A(n—r—1)1 ^ 

A.Au a O 


(3.28) 


Here Am = Am + um = um + '&wVm and Tmn, ^mn are their field strengths. C 2 (r) 

is a constant. One can check that no other terms are allowed. Note that since entropy current 
is not a direct observable, only its divergence is, we are allowed to include gauge non-invariant 
terms in as long as the divergence is gauge invariant. Computing divergence of eqn. (3.28) 
we can obtain. 




n—1 


+Er 7 )'"-> 

r=l 


—V M<^s{r) — 


E + P — dwR 


R 


C2,{r) - ^T^C2,{r+l) ) 


Q, 


+ ( - C'2,(r+1) ) {^MNU^ “ 


(3.29) 


On the other hand divergence of canonical entropy current through eqns. (3.12) and (3.14) is 
given as, 

n—1 




r=0 


A 


(3.30) 


Combining the two pieces and demanding VmJs ^ 0, we find a consistency condition in 
entropy current that a)s(r) = d)s(r)('d) must not be a function of w^v. Prom here parity odd 
transport coefficients in eqn. (3.13) are determined to be. 


We(0) = (i?C'i^(0) - ^vC2^^^y)) , £l;q(0) = 


(3.31) 


and for r / 0, 


^e{r) = + 


E + P — -dwR 
R 


C'2,(r) - '&^C‘ 


a.C-hi) 


Q. 


^q{r) — ( ^C'2^(r) — C'2^(r+1) 


(3.32) 


^The C’s in this expression are arbitrary, and a priori have no connection to the C’s introduced in previous 
sections. However, as is suggested by the notation, both will eventually turn out to be the same quantities in 
constitutive relations. Also, the entropy current need not be gauge invariant. 
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where = w‘^s{r){'&) is an arbitrary function of 'd. Compared to equilibrium partition 

function constraints eqn. (3.19), we have one additional constraint, 

(3.33) 

and one less constraint: Ci is not a constant but a function of "d. In appendix (A) we 
show that on introducing torsion, entropy current positivity will indeed set C*! to be a 
constant. 


<^2,(0) — 0) 


Leading Order Parity Odd Sector (for d = 2n): Now we perform a similar analysis for 
even dimensional parity odd sector. Similar to odd dimensional case, here also we can have 
terms in whose divergence does not have any pure derivative terms^^. 




(r) 


r=0 


uAT An^ ^ 


n—1 


+E("")'(T 


r=0 


^s{r)^N^ • (3.34) 


One can check that any other term if included, will give pure derivative terms in divergence. 
Divergence of this object can be computed fairly easily to be. 


VMTf = 0^ (”)/(. 


r) 


r=0 


^dP d IdP d IdP d 
^ 'ddR'^ ^ 'ddQ'^ 


>2,(r) 


n—1 


I /n-l\iMN 

\ r rw 


r=0 


^ X ^ r. , + P -dwR d ^ ^ c \ O ,00 _ 

K-w{r) nv ^^S2^[r+l) ] Qu'^OnW 


dd 


dw 


dd' 


4“ 1 ^su(r) ^s{r) 4“ ^ 


E + P-dwR d d 

S 2 ,{r) - ”^^‘ 52 ,(r+l) I OMdOnV 


dR 


du 


d 


d 


+ dM^d^w + n {^^dMS 2 ^(r) - dMS 2 ^{r+i)^ {Enru^ “ ddN^) 

(3.35) 

On the other hand divergence of canonical entropy current through eqns. (3.12) and (3.15) 
is. 


n—1 


- = -^Y.OA)hr) + E (" 4 )' 


-l^ iMN 
(r) 


r=0 


r=0 


Kmir)^dMddNW + Kuir)^dMddNl^ 


'd 


n—1 


+E (”4)'r 


(r) 


r=0 


4“ 4“ id^q{T^ 4“ ^efr)) ^^Ivld 


{Enru — ddNi^) 


(3.36) 


^^The 5’s in this expression are arbitrary, and a priori have no connection to the S’s introduced in previous 
sections. However, as is suggested by the notation, both will eventually turn out to be the same quantities in 
constitutive relations. 
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Combining the two pieces and demanding S/mJ^ > 0, we get a consistency condition on 
entropy current, 

(3.37) 

whose most generic solution is, 

d d d 

^svj{r,s) /l ; ^su{r,s) Qp ’ (3.38) 

for some functions /i(d, ro, z^), f 2 {'d^u). We define, 

^s{r) “1“ /l■ (3.39) 

Expressed in these variables, one can check that entropy current positivity gives all the par¬ 
tition function constraints eqn. (3.24), except the expression for modifies to, 


^eu(r) 




dud^ 


/ 2 ('d, v) + "d—-|- 


E + P- dwR d 


d 


R 


^‘52,(,)-z9z.^52,(,+i) 


, (3.40) 


and in addition we get, 


(3.41) 


Like even dimensional case, we again see that we get an additional constraint through en¬ 
tropy current, but one constraint turns out to be weaker. Equilibrium partition function 
a^a^ /2(^; ^) = 0 which entropy current fails to do. In appendix (A) we will show that 
introducing torsion remedies this situation. 


^g(r) 4“ 


Qd d 


Rdi9 


3.4 Recap 

In this section we summarize the results for leading derivative order null fluid in mass frame, 
taking into account constraints from equilibrium partition function and second law of ther¬ 
modynamics. The constitutive relations for null fluid are given in terms of fluid variables z?, 
tu, V, , 

J,MN ^ + ppMN -h 

= Qu^ + T^, (3.42) 

where P, R, E, Q are thermodynamic pressure, mass density, energy density and charge den¬ 
sity expressed as functions of ??, w, v. These constitutive relations follow the conservation 
laws, 

VmJ^^ = 0. (3.43) 
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In odd number of dimensions (d = 2n — 1), the form of hydrodynamic frame invariant correc¬ 
tions to leading order in derivatives are given as, 

n—1 

= KeP^^dN'd + 'dKqP^^ {PnRU^ “ d^Arl/) ^ , 


r=0 

n—1 


= KqP^^dN^ + C7qP^^ {PnRU^ " ^q{r)ltry ( 3 - 44 ) 


r=0 


where transport coefficients r] (shear viscosity), C (bulk viscosity), aq (electric conductivity) 
are some non-negative, (thermal conductivity) is a non-positive and Kg (thermo-electric 
coefficient) is an arbitrary function of zu, v. Parity-odd transport coefficients (Hall con¬ 
ductivities) are however completely determined upto some constants as, 


E + P — -dwR 

R 


d’e(r) = + 

^q(r) ^2,(r’-|-l)^ • 


C'2,(r) - '&^C‘ 


a.P-hi) 


(3.45) 


where C’s are some arbitrary constants, and ^ 2 ,( 0 ) = 0- even number of dimensions 
[d = 2n) however the corrections are given as. 


n—1 


Y^mn ^ _ 


ce + E OcV,i,„ , 

r=0 \ r=0 / 


‘r=0 


‘r=0 


n—1 




r=0 


n—1 


n—1 


r=0 


r=0 


pM = ( P^^Kq + E I + ( aqP^^ + E C~r"Kf^<lir)j {J^NRU^ " 

(3.46) 


n—1 


-h dn E Tr ( -^dNS2,{r) - dN‘S2,{r+l) 

r=0 ^ 

where we have made following redefinitions with respect to eqn. (3.15), 

^e(r) ^ ^e(r) + 


^ ^ E + P- dwR d ^ a ^ c 

+ p - du—S2^^r+l) 


^q{r) ^ ^q{r) 4“ 


Q d 


d 




(3.47) 


The transport coefficients in parity even sector are same as before; however parity-odd trans¬ 
port coefficients i)(r) (Hall viscosity), (thermal Hall conductivity), Kqi^r) (thermo-electric 
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Hall coefficient), ^q[r) (electric Hall conductivity), and S 2 ^{r) are some arbitrary func¬ 

tions of 'd, w, u. Finally C(r) is determined as. 


C(r) 


& 


.dP d dP d dP d 






dEdd dRdw dQ du 


^2,M- 


(3.48) 


All the constitutive relations satisfy the physical requirements of existence of an equilibrium 
partition function and entropy current. To leading order in derivatives they are given as, in 
odd number of dimensions (d = 2n — 1), 


n—1 


r=0 

n—1 

(L„) + E * [{"{”;‘)Cl.(,)dV - (4i)Cj,|,+i,.4} A « A A f! 


jM ^ jM 


* Ar AA(n—r—1)1 ^ 


r=0 


(3.49) 


and in even number of dimensions (d = 2n), 


/ .. n—1 n 

{dA‘} -EO'oM'S-dr) 

. ° r=0 r=0 


jM _ jM I 

s{can) ' ^ 


n—1 


n 


E (E)‘5i.M V A di9 + {DS2M^ 


f.Ar AA(n—r—1) 


r=0 


r=0 


M 


(3.50) 


While writing entropy current from eqn. (3.34), some total derivative terms have been dropped, 
as they will have zero divergence. We have included the additional constraints coming from 
entropy current analysis while writing the partition function and vice versa. This finishes the 
discussion of null fluid upto leading order in derivatives in arbitrary number of dimensions. 
Next we turn on to study the light cone reduction and how to get Galilean fluids via reduction 
of a null fluid. 


4 Light Cone Reduction 

We want to study a Galilean system in (d -|- 1) dimensions. So we essentially want to com- 
pactify the V direction as M.(^d+ 2 ) = 5'^ x But, V is null, and thus is transverse to 

itself, so it is not possible make such decomposition uniquely. It is therefore convenient to 
introduce another vector field (we call time field) T = T^du, which can be used to define a 
unique decomposition Ai(^d+ 2 ) = *5)/ x x where ^(d) is the set of vectors transverse 
to both V and T. The time field T provides a reference frame for compactified Galilean 
theory. We formally define light cone reduction as this choice of frame and subsequent com- 
pactification. 
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Using T we can define another null field orthonormal to V, 


^ 

(r) t^Vn 

which satisfies = 0; V(x)^M = — 


1, and a projection operator transverse to V and 
-^fd) = V-TV : € Af (,+2)} • (4.2) 


Since the choice of T is arbitrary and does not have any physical significance, null theories 
are invariant under an arbitrary redefinition of T —)• T', which we parametrize as, 

tM = a [T^ - T^Vac!/’^] , (4.3) 

where = 0 and a G M. One can check that inverse transformation is simply a —>• 

1/a, and —>■ Our funny parametrization has a benefit that under T redefinition, 

transformation of and only depends on 

^ + (4.4) 

where So our light cone reduced theory is described on a compactified 

null background A4(^d+2) with isometry ipv = : ^{V)}: and a time field T^, modded by 

diffeomorphisms, gauge transformation and T redefinition. In light cone reduction approach 
however, we need not worry too much about T redefinition. Since the original theory on 
null background did not depend on T, so the reduced theory will also be invariant under its 
redefinition automatically. 


4.1 Newton-Cartan Backgrounds by Light Cone Reduction 

It is easy to see how the Newton-Cartan structure comes out by light cone reduction. We 
identify x as the degenerate Newton-Cartan (NC) manifold. Without 

loss of generality we can choose a basis in the original manifold A4(rf+2) such 

that ipv = {y = = 0}. will then provide a basis on NC manifold This 

mechanism to generate NC manifold via null reduction was first found in [12] and has been 
further developed in [6, 13]. 

Reduction of Background Fields: We can decompose background fields according to 
this choice of basis as, 



- 25 - 




P{T)MN 


(4.6) 


/^O o\ MTV ^ 

1^0 p^,) ’ 1 ^ p^^'^B, ) 

such that 

n^v^ = 1, = 0, n^p^" = 0, P^ipP^'' + n^v"' = 6^"'. (4.7) 

This is the well known Newton-Cartan structure. The T-redefinition invariance introduced 
above becomes Milne boosts on NC manifold, 

+ + - 2n(^Vi^)+n^ni.Vi^. (4.8) 

from where we can identify as NC frame velocity. Similarly the Levi-Civita connection 
can be reduced in NC basis as, 

^ pu ~ ^ T 2^ ^ i^pPpi' T duPpp — dpP^y) + ^{T)a{v^v)P 1 

^~pu = P\{p^u)V^ - (4.9) 

and all other components zero. Here we have identified as Newton-Cartan connection 
and denoted respective covariant derivative by V^. In eqn. (4.9) we have treated derivative 
of frame velocity ^pv'^ as an independent variable, and have defined frame vorticity in terms 
of it as ^(T)pu = “^Palu^p]V^■ The covariant derivative can be checked to act on NC structure 
appropriately, 

'^pny = 0, '^ppf'' = 0, '^pPvp = -2n(^pp)^V^u°'. (4.10) 

We define the spatial volume element on NC manifold as, 

^pu... ^ ^MNpy...y^y^^^^ ^ -e-P^^^-Up, Spy... = PppPy„ . . . eP^-. (4.11) 

and * as the Hodge duality operation associated with it. The notations and conventions on 
differential forms can be found in appendix (D). Finally the only surviving components of 
gauge field strength are Tpy, which can be identified as NC gauge field strength. We can 
further decompose Fpy and ^(T)pv into, 

^py = 2e[^nj,] + Ppy, ^{T)py = —‘^0;{T)[p'P‘'y] + ^(T)py (4-12) 

All the introduced tensors are transverse to v^. Here is the electric field while /3py is the 
dual magnetic field defined with respect to frame T. Similarly a(T)/i is frame acceleration and 
^{T)py is spatial frame vorticity. We can similarly define the fluid acceleration and vorticity 
as well, which will be used later in § 5. 
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Reduction of Currents: We decompose the currents of null theory on M.{^d+ 2 ) — 5'y x 
X as, 

+ 2etotV^^V^] + 24 ^V(Jj + + OiV^V^, 

+ e2V^. (4.13) 

All the introduced tensors are projected along , Note that redefinitions eqn. (2.8) can 
be used to get rid of 6 terms. Choosing Newton-Cartan basis we identify, as stress-energy 
tensor, p,jp as mass density and current, etot, je as energy density and current, and q,jq as 
charge density and current, as seen by frame T. Under a finite T redehnition eqn. (4.3) they 
transform as, 

p ^ p, fp ^ fp - pr, ^ + prr, 

etot etot - + \p'4’‘^: Je {je “ (^totV) “ {t^'' “ Vjp) 'iju + {jp “ , 

q^q, j^ j^ - qi)^- (4.14) 

Note that we have used the same time field (reference frame) to decompose currents/densities 
as well as background fields. It is sometimes required to define background in one reference 
frame (e.g. lab frame) but currents and densities in some other reference frame (e.g. co¬ 
moving frame). One can merely perform a T redefinition on various quantities noted above 
and gain the desired result. 

Reduction of Ward Identities: On the decomposition A4(£^_|_2) = Sy x x 
background field content is Vm, V(t)Mj P(t)mn and Am-, so any physical theory should 
be described by a partition function lU[yM) .P{T)MAf) Using the current redefini¬ 

tions eqn. (2.7) for null backgrounds, we can parametrize the variation of partition function 
as, 


6W= f {dx^} [{etotVff)+j^ 


SVm + 
1 


\P^(T)+Jp 


6V( 


(T)M 


I I I \ CD I 

+ ( P UT) + ) oP{t)mN + 


\qV{T) + 3q 


5Am 


(4.15) 


The same partition function should also be gained by directly reducing the relativistic parti¬ 
tion function eqn. (2.4). This will render the quantities in above partition function to be same 
as the ones defined in eqn. (4.13), and in addition . The latter constraint is Ward 

identity of T redefinition, i.e. can be gained by demanding partition function eqn. (4.15) to 
be invariant under T redefinition. 
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Choosing NC basis the partition function variation eqn. (4.15) can be decomposed to, 


5W = f {dx^} ^det(p^,^ + n^Uu) 


- (etotv^ + if) + {pv^^ + if) 6Bfj 


+ P^V' + 5p^,y + {qv^ + if) 5A^ 


(4.16) 


Symmetry data of the light cone reduced theory is, 


( 4 . 17 ) 

where we identify as mass parameter, as time translation parameter, as space 
translation parameter, as gauge parameter, and as Milne boost parameter. The 
respective Ward identities can be found using eqn. (4.16) or directly reducing the null Ward 
identities eqn. (2.5), 


Mass Conservation: 
Energy Conservation: 
Momentum Conservation: 
Charge Conservation: 
Milne Identity: 


Vm {pv^ + jji)=o, 

(etotv^+j^) = - {v^py + t^y) 

(u^p- + t^-) = [qe'^ + /3yf] - (pu^ + jf) 

(c?^;^+if) =0, 

fp=P^- (4.18) 


First terms in the RHS of energy and momentum conservation equations are work done 
and Lorentz force due to electromagnetic fields, while the last terms are pseudo-energy and 
pseudo-force due to spacetime dependence of the frame velocity. As we already mentioned, 
the Milne identity is trivial in theories obtained by reduction. 


All these results here have been mentioned in Newton-Cartan notation, which is a nice covari¬ 
ant formalism for Galilean physics. However for familiarity and to build intuition, we have 
given all these results in conventional non-covariant notation as well in appendix (B). 


4.2 Equilibrium on Newton-Cartan Backgrounds 


From the perspective of Galilean theories, equilibrium is defined by a preferred reference frame 
(or time field) K with respect to which system does not evolve in time. This can be achieved 
by reducing null theories at equilibrium, and identify the timelike isometry ■i/’A' with preferred 
reference frame in the Galilean theory. Hence the variation of eqb. partition function in local 
rest of reference frame K is essentially same as the null fluid eqn. (2.14) written in terms of 
Galilean quantities. 


= / {dx*} 


1 


1 


+ ( -—foe + ^ofop + T^ofoq ) — 


1 






A. 


+ ( qSl^o + '^foq^^i 


(4.19) 
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and hence, 




Here all the observables are defined as seen by reference frame K, and are denoted by a 
subscript o- These will reduce to the expected relations eqn. (1.3) in fiat space, i.e., 


1^0 — 1 , Qij — ^ij'! '^o — — -^i — — 0 * 


(4.21) 


Since we have fixed that T-redefinition symmetry by choosing a preferred reference frame K, 
the corresponding EOM does not show up. Consequently momentum current does not 
appear in the partition function, and can be found by using the missed EOM. In equilibrium 
configuration, null fluid and Galilean fluid have same field content and symmetries, so we 
expect the eqb. partition function to also be the same. To ideal order eqn. (2.21) it will 
identify p,e,q with thermodynamic functions R,E,Q, and hence will give physical interpre¬ 
tation to thermodynamics of null theories in terms of Galilean physics. In hydrodynamic 
description, at further derivative orders also, it will give physical interpretation to various 
transport coefficients and constraints of null fluid. 

5 Galilean Hydrodynamics 

Having discussed the light cone reduction of generic null theories in last section, we can 
straight away perform light cone reduction of null fluids in § 3 and hope to get Galilean 
fluids. If we look at eqn. (3.1) closely, it is already nicely organized in local rest frame of the 
fluid (defined by = u^). One just needs to apply a T-redefinition to it with = —u^ 
to get densities and currents in a generic reference frame. 



where we have identified 


p = TZ, q= Q, e = £, 


<;-M = F, = (T - P)p^'' + (5.2) 

and = u^u^. Similarly entropy current of the Galilean fluid can be found to be. 











which follows second law of thermodynamics 


(5-4) 

Choosing ‘mass frame’ in the null fluid, which is the most natural frame from a non-relativistic 
perspective, will switch off and hence Galilean mass current will not undergo any 

dissipation. The identifications for mass frame are given by eqn. (5.2) can be read out in 
terms of frame invariants as, 

p = R, q = Q, e = E, s = S-Ts_, u^ = u'^pj^, 

= 0, = T^, = S^, (5.5) 

and in turn constitutive relations become, 

=p^^ = Ruf^, = Ruf^u’^ + Pp^^'^ + = Qu^ + c;^, 

^tot = E —RE^, = ^E + P + E E TT^’^Ui,, 

s = S- T,_, = Su>^ E - Ts-nE (5.6) 

These are the standard Galilean constitutive relations, written in Newton-Cartan basis. We 
will present all these expressions in conventional non-covariant basis in appendix (B) for the 
benefit of readers not comfortable with Newton-Cartan formalism. 

Having obtained the general picture, we can now deduce constitutive relations for a Galilean 
fluid upto leading order in derivative expansion, using the corresponding null fluid results in 
§ 3.4. In table (2) we have mentioned light cone reduction of all the leading order data to 
get Newton-Cartan data. Having done so, rest of the algebra is essentially trivial. In the 
following we will work in mass frame exclusively. 

Even Dimensional Galilean Fluids: Using reduction of data enlisted in table (2), we 
can read out the even dimensional {d = 2n — 1) constitutive relations from § 3.4, 

n—1 

= l^qP^'^dud -h (JqP^^^ {RypU^ - ddyU) -b ^ (^~^) Coq{r)l(ry 

r=0 

n—1 

= KeP^'^dyd E dKqjE"' {TvpVE - ddyV) E ^ 

r=D 

where transport coefficients p (shear viscosity), C (bulk viscosity), aq (electric conductivity) 
are some non-negative, (thermal conductivity) is a non-positive and Kq (thermo-electric 
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Table 2: Leading Derivative Order Data for Galilean Flnid 


Null Fluid Data 

Newton-Cartan Data 

Parity Even 

0 = Vmu^ 

0 = ^puP 

P^^dnd, P^^dNW, P^^dNV 

pP'-'dyd, pP'^diyZu, pP'^d^y 

P^^ {Pnru^ - ddny) 

pP^ (P^pup - dd^v) 

aMN = 2 P^^P^^V^rUs) - IP^^Q 

III 

q 

1 

a-|bo 

0 

Parity Odd — Odd Dimensions (d = 2n — 1) 

iM 

hiL 

” * (P A u A A ^ 

r=0 L J 


Parity Odd — Even Dimensions (d = 2n) 

hr) 

” =*rPA-uA:r^’'A 

r=0 L J 

hr) 



^(Z)^ud, 

lf^{^ {Pnru^ - ddNR) 

ip;^^ {p,pu- - dd,v) 

,R(M N) 

hr) ^R 

f{R u)(T 
‘‘(r) PP dpa 

where, 

=*[v A 

ki 

= * [/3^^ 


coefficient) is an arbitrary function of i?, w, v. Parity-odd transport coefficients (Hall con¬ 
ductivities) are however completely determined upto some constants as, 


E + P — -dwR 

R 


^e{r) = + 

^q{r) ~ , 


^2,{r) I ) 


(5.8) 


where C’s are some arbitrary constants, and ^^ 2 ,( 0 ) = 0. As a special case one can obtain the 
4 dimensional {d = 3) results which will only affect the parity odd sector, 

+ aqp^'' {PupUp - dd^v) + 

= KepP'^dud + dKqpP" {PypUP - ddyV) + oje{0)UJ^ + , (5.9) 

where, 

= ^io) = BP = = ^e'^PP^p, (5.10) 
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are vorticity and gauge magnetic fields respectively, and, 

^q{0) = —2??C'2,(i), £Uq(i) = 2?? — <^2,(2)^ • 


E + P — dwR 
R 


^2,(1) “ "^^^2^(2) I , 


(5.11) 


Odd Dimensional Galilean Fluids: Using reduction of data enlisted in table (2), we can 
read out the odd dimensional {d = 2n) constitutive relations from § 3.4, 


n—1 




(r) , 


r=0 


r=0 


< = Tr + [P^'^^q + Y ("r ^ W^<?(U j 

V r=0 / V r=0 / 

PdnY rYlr) {%dyS2,(r) “ , 

r=0 ^ ^ 

n—1 \ / n—1 \ 

+ E iYKY(r) ^ - Y {YK)^(r) 

r=0 J V r=0 / 


= P^^P 


n—1 


+ ^^E (E)^M UduSi,(^r) + ^ - ^iyd,S2,(r+l) ) • (5.12) 

r=0 


The transport coefficients in parity even sector are same as before; however parity-odd trans¬ 
port coefficients (Hall viscosity), (thermal Hall conductivity), K,q(r) (thermo-electric 
Hall coefficient), cT^(p) (electric Hall conductivity), 5^ (p) and S 2 ^{r) are some arbitrary func¬ 
tions of T?, w, V. Finally ("(p) is determined as, 


C(r) = 


■ gSP d dP d dP d 
dE d'd dR dvo dQ du 




(5.13) 


As a special case we would like to write down the 3 dimensional results, 

(C0 + C^o; + (bB) , 

^q = {p^^'f^q + Su'd + {P^'"CPq + e^'^dg) {PypU^ - 'OdyV) 

+ (^^^dyS2^{Q) — 'dduS2^(i^ , 

= {p^’^Ke + e^'^Ke) dyd + d {p^'^Kq - Kq) {PypU^ - ddyV) 


P 'de^^ (Q) -|- 


F -|- P — dwR 


duS2^{o) - dvdyS 


P‘P2,(l) I ) 


where. 


^ — ^(0) — 2^^ 


B = hi) = 


(5.14) 

(5.15) 


- 32 - 








are again the vorticity and gauge magnetic fields and we have renamed dg = dg(o)) = ^e(o)) 
Kg = Kg(o)! 1? = f?(o)i Co; = C(o)i Cb = C(i)- The 3 dimensional Galilean fluid was also studied 
by [ 2 ], however we find certain discrepancies in their and our results. A detailed comparison 
has been provided in appendix (C). 

Before closing this discussion we would like to note that [7] also constructed an equilibrium 
partition function and entropy current for an uncharged 3 and 4 dimensional Galilean fluid, 
and used it to constraint the respective constitutive relations. By switching off the charge 
sector and setting <7 = 2,3 we see that we trivially recover their results. 

This finishes our discussion of (non-anomalous) constitutive relations of a Galilean fluid upto 
leading order in derivatives in arbitrary number of dimensions, obtained by light cone reduc¬ 
tion of a null fluid. Unlike the hydrodynamic reductions before this work [14, 15], there is no 
non-trivial mapping between the relativistic (null) fluid and the Galilean fluid. In fact term 
by term, null fluid constitutive relations are same as Galilean constitutive relations. Same is 
true for thermodynamics, entropy current and the equilibrium partition function as well. We 
deduce that we can see null fluid as Galilean fluid written in extended space representation. 
Many aspects of it are already hinted by extended space construction of [2]. In next section 
we extend this approach to study effect of 17(1) anomaly on fluid transport. 

6 Anomalies 

Upto this point we have studied hydrodynamics on non-anomalous null/Galilean backgrounds. 
In this section we want to explore if the null background construction can also be used 
to introduce 17(1) anomaly in Galilean theories^^. Later we will hnd how constraints of 
Galilean fluid modify in presence of anomalies. We use the anomaly inflow mechanism of 
usual relativistic backgrounds to achieve this goal, with appropriate modifications due to the 
null structure of the background. 

Gonsider a bulk manifold ^(^+ 3 ), on whose boundary A 4 (rf+ 2 ) our theory of interest, i.e. 
null fluid lives. Indices on are denoted with a bar M,N.... We define ^(^+ 3 ) also 

as a null background, with respective fields G'Miv ^ compatible null isometry^^ 

ipv = {17 = d]^, A(y^ = 0 }, such that transverse components of all these fields vanish at 

boundary. We can define respective fields on M.(^d+ 2 ) by pulling back the bulk helds, which 
gives it a null background structure. 

Galilean anomalies consi<lered in [18] are different than what we are considering in this paper, because our 
background field content does not match that of [18] after reduction (we have chosen A- =0). A detailed 
comparison of these issues along with an extension to non-abelian and gravitational anomalies will shortly 
appear in a companion paper [19]. 

^"^We would like to mention that this construction only seems to work when we set A(y) = 0. We give more 
reasoning in this regard in a companion paper. 
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We start with the assumption that full theory on ^(^+ 3 ) U Wl(d+ 2 ) described by a partition 
function W is gauge invariant. Most generic such partition function can be decomposed into 
a bulk and a boundary piece, 

W = W[M^d^2)] + ( 6 . 1 ) 

which individually are not gauge invariant. Here W is the partition function of the boundary 
null theory which is anomalous, i.e. is not gauge invariant. Whuik on the other hand is a 
pure bulk piece whose gauge variation must be a boundary term. While constructing Whuik 
out of background fields, we can let go of any terms which are gauge invariant upto a total 
derivative (we can always redefine W to absorb this total derivative term at the boundary), 
as they will not induce any anomalies in the boundary theory. Hence allowed Whuik can be 
written as integration of a full rank form, 

Whuik = [ ( 6 . 2 ) 

such that has an exact gauge variation = dG(^), and it must not be symme¬ 
try invariant upto an exact form^^. In usual relativistic theories, can only be written 

in odd bulk dimensions (d = 2 n — 2 ), and is given by the Chern-Simons form 

/g+d = (6.3) 

However for null backgrounds, this term identically vanishes, as it is a full rank form but 
does not have any component along V. We are therefore forced to modify by adding 

some term which has non-vanishing component along V. We do it by choosing an arbitrary 
time-field T and use it to define a conjugate null field H(t). Now we can define an analogue 
of Chern-Simons form, but in even bulk dimensions (d = 2n — 1), 

j{2n+2) ^ ^ A ^ . (6.4) 

We need to check if it fits our requirements. We will leave it for the readers to convince 
themselves that this expression cannot be transformed into a symmetry invariant term by 
adding an exact form. For the other criteria we need to compute its gauge variation, 

(6.5) 

Last term vanishes as it has again no component along V, thus we verify that gauge variation 
of /( 2 ^+ 2 ) is a boundary term. It is important to note that while we have used V(j'^ to define 
j(2n+2)^ it is invariant under T-redefinition. One can check there does not exist any other 
term which meets these criteria. Hence contrary to usual relativistic backgrounds, here we 
can only define /( 2 ^+ 2 ) even bulk dimensions. It means that only odd dimensional null 
backgrounds (the one at the boundary) and hence even dimensional Galilean backgrounds 
(that we get by reducing the boundary null theory) can be anomalous, which is what we 
expect. 

/ 7^ /'-t- dJf, for some gauge invariant 
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6.1 Anomalous Ward Identities 


In this section, we present the modified Ward identities in presence of U{1) anomaly. In 
presence of anomalies, variation of boundary partition function W generates consistent cur¬ 
rents and J^ns which are not gauge invariant. Varying /(2"+2) write down 

variation of the full partition function, 


6W = 


j {n+l)C^‘^'^hA^V^T)^^^'"+ [ {dx^} VG 


1 . 




M 1 


, ( 6 . 6 ) 


where we have defined the covariant currents, 


nMAf 


= T, 


:MN 


* \V(r) A ^ A 


1 M 


(6.7) 


Since the full partition function is gauge invariant, and we see that the bulk piece is manifestly 
gauge invariant, therefore covariant currents must also be gauge invariant. Demanding W to 
be gauge invariant we can get the Ward identities in the boundary theory, 

= -{n + * [F(t) A . (6.8) 


We observe that the system exhibits U{1) anomaly. 


6.2 Anomalous Equilibrium Partition Function 


In our earlier discussion on equilibrium in § 2.2, we wrote the most generic equilibrium 
partition function as a gauge invariant scalar. Now we need to modify this partition function 
appropriately with a gauge non-invariant piece to account for anomaly. We decompose the 
equilibrium partition function into, 


= W!,ts + W^tm- (6.9) 

Here WcSns is the most generic gauge invariant partition function which can be written out 
of background helds, which has been discussed thoroughly in preceding sections. Wanom on 
the other hand is completely determined in terms of anomaly coefficient We suggest 

its explicit form to be. 


JM^2n + l) \ ^ J 


( 6 . 10 ) 

which generates correct anomalies. Recall that we are allowed to use any arbitrary time 
field T to specify anomalies; in equilibrium configuration T = K is the most natural choice. 
Performing on shell variation (dV = 0) of Wa^am, 


^ ' nnr^m 


IM 


-n(n + ^V ^ V^k) A 


(2n+l) 


IM 


(2n+l) 


{(n + l)??oJ/oF A V^K) A - V^k) A ^ A , (6.11) 
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and using eqn. (6.7), we can find find the covariant anomalous currents at equilibrium, 


J^anom = n{n + * 

T^anon. = n{n+l)&^'-^lvl^ 


V A V(K) A 


M 


1 


F A y (K) A 


( 6 . 12 ) 


One can check that these currents identically satisfy the anomalous conservation equations. 
Note that these currents are also to be supplemented with the non-anomalous pieces discussed 
in previous sections. In the local rest of reference frame K, the equilibrium partition function 
can be expressed in Kaluza-Klein notation, 

W2tn = [ A f(dA)''(’^-^) + \{n- l)uoMa A (dA)^^’^-^)) . (6.13) 

JM^2n-l) V ^ / 

Above we have left more than one powers of da, as they do not contribute on-shell. Varying 
it we can find the anomalous contribution to Galilean currents; only non-trivial contributions 
are given by. 


fq,anom = n{n + l)C^‘^'^'^'doVo* [(dA)^^’" 




(6.14) 

When generating constitutive relations of a Galilean fluid using equilibrium partition function, 
above results are naturally written in equilibrium hydrodynamic frame. Interestingly mass 
current does not get any anomalous correction, hence these results are automatically written 
in mass frame as well. Correspondingly the constraints of parity-odd sector in odd spatial 
dimensions (d = 2n — 1) eqn. (3.19) modify to include contribution from anomalies, 


^q(n-l) = — C2^(n) + {n + l)zzC^^"'^^ . (6.15) 


Note that only r = n — 1 component of eqn. (3.19) is modified, while other constraints remain 
unchanged. 


6.3 Anomalous Entropy Cnrrent 

In this section, we shall try to get the anomalous contribution to constitutive relations found 
in last sub-section, using second law constraint. In presence of anomaly, the canonical entropy 
current divergence eqn. (3.12) will get modified to, 

= -n^^VMUiv - {Tmnu^ - 

+ i?i/(n + l)C(2 ’")*[F(T) a:F^”] . (6.16) 
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Using leading order parity-odd {d = 2n — 1) constitutive relations eqn. (3.14) we can evaluate 
it to get, 


{Tmnu^ — "ddu^) 

- ^vn{n + ■ (6.17) 

Clearly it will modify the constraints eqn. (6.15) only for r = n — 1. Plugging in the expression 
for V from section § 3.3.2, we will reproduce expression for in eqn. (6.15), and 

get the following differential equations for 


iVK) 


-^Cbs(n-i) = 0, -^Cbsin-i) = n{n + (6.18) 

Last equation will imply, 

u^sin-i) = ^n{n+l)diy'^C^‘^^^+f{d), -^Cjs(n-i) = ^n{n+l)iy‘^C^‘^'^'^+Ci^(^n-i)i'&)- (6-19) 

This gives the remaining Co^(^n-i) constraint in eqn. (6.15), except that Ci (^n-i) is an arbitrary 
function of "d similar to what we saw in § 3.3.2. This can be remedied by putting in torsion 
to relax T-Lmn = 0 condition, as we shall present in the appendix (A). 


7 Discussion 

In this work we have proposed an innovative and interesting approach to construct the equilib¬ 
rium partition function and constitutive relations of a Galilean fluid in arbitrary dimensions, 
starting from a relativistic system, namely null fluid. The basic idea of this construction has 
already been presented in a previous paper [1]; here we have generalized it to include a global 
(anomalous) U{1) current. The beauty and importance of our approach lies in the construc¬ 
tion of null fluids. We have showed that the symmetries and background field content of a 
Galilean theory is exactly captured by a theory defined on null background. This motivates 
us to dehne a theory of hydrodynamics on null backgrounds (i.e. null fluid) from scratch, and 
use it to derive constitutive relations of a Galilean fluid. 

Although, the main aim of this construction has been to write down leading order constitutive 
relations of a Galilean fluid (in presence of anomaly), but in the process we learnt that null 
fluids can be considered as a robust stage to study properties of the most generic Galilean fluid. 
We found an exact one to one correspondence (not corrected order by order in derivatives) 
between all aspects of a null fluid and a Galilean fluid, but more than that, the actual map of 
this correspondence is essentially trivial. Our approach has been to study the null fluid itself 
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as an independent theory, and later exploit the triviality of this map to say something useful 
about the Galilean fluids. 

The triviality of this map is an important feature. Works in past in this direction, including 
our own in [15], have at best found a mapping between usual relativistic fluids and (a subset 
of) Galilean fluids, which has to be corrected order by order in the derivative expansion 
(not to forget the faulty thermodynamics it endows to the Galilean fluid), and not much 
useful could be said thereof. Note that in the current work however, our mapping is exact 
to all orders in the derivative expansion, which enables us to directly use much sophisticated 
relativistic machinery to study non-relativistic fluids and hence is very interesting in realistic 
scenarios. 

We also found that null backgrounds allow us to introduce U{\) anomalies in an odd¬ 
dimensional null theory (i.e. an even dimensional Galilean theory) and forbid them in an 
even dimensional null theory (i.e. an odd dimensional Galilean theory). This is in sharp 
contrast with the usual relativistic results, where it is the other way round. However from 
the perspective of Galilean fluids it is pretty natural. A generalization of this construction 
to include non-abelian and gravitational anomalies will be presented in a companion paper 
[19]. 

Apart from these anomalous terms, we also have other parity-odd terms in the constitu¬ 
tive relations, for both, even and odd dimensional Galilean fluid. The study of parity-odd 
non-relativistic hydrodynamics has become a very fascinating topic in recent years. Fluid 
consisting of chiral molecules breaks parity at the microscopic level. This kind of fluid plays 
an important role in many biochemical processes, where only the molecule with right chirality 
can fit into a protein. Therefore to model such a fluid, we would be forced to add parity-odd 
terms in the constitutive relations. Our construction gives a consistent way to obtain the 
possible parity-odd terms in a Galilean fluid, at any desired order in the derivative expansion. 
It would be very interesting and important to understand the effect of these terms in some 
practical examples. 

Finally, we would like to make some comments on physical aspects of the null fluid. Although 
we construct the null fluid dynamics and show that it is in one to one correspondence to a lower 
dimensional Galilean fluid; the physical significance of this null fluid itself is not yet clear. 
Presence of an extra background field, i.e. the null Killing vector, has allowed us to introduce 
a set of new transport coefficients (e.g. R) in the null fluid constitutive relations, as compared 
to a usual relativistic fluid. The physical meaning of these new transport coefficients, which 
come coupled to the null vector, becomes clear only once we identify the dynamics of null 
fluid with that of a Galilean fluid living in one lower dimension. Hence, at this stage, the 
correct physical interpretation of the null fluid appears to be that it is a particular embedding 
of the Galilean fluid in a spacetime with one higher dimension. This approach is more in lines 
with the axiomatic approach of defining a Galilean fluid, but has the benefit that we have all 
the well-developed machinery of relativistic physics at our disposal. 
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A Minimal Torsion Model 

In the main text we alluded at more that one instances that entropy current fails to cap¬ 
ture all the constraints a null/Galilean fluid should follow on torsionless backgrounds. In 
this appendix we provide heuristic arguments for consistency between entropy current and 
equilibrium partition function constraints in parity-odd sector by introducing ‘minimal tor¬ 
sion’. We introduce a torsion in the connection, = — \u^'Hmn-, so the full connection 

becomes, 

k'Vs' = 2*^^^ {Qm^ms + QsGnm — QnGms) — 2 — umT-Ls^ — usT-LjJ^) ■ (A.l) 

One can check that with this connection, VmV^ = 0 does not require T-Lmn to be zero. In 
fact, for our purposes it suffices to choose H. = dl^ = V/\X where is a projected vector 
such that dA = 0. We call this minimal torsion. For spinless theories, constitutive relations 
remain the same except Vm —^ ^m — Xm- We can define an entropy current as before, 

whose canonical part will now have divergence, 

^^MJ^can) = (am - {XmNU^ " • (A.2) 

We will now write the most generic constitutive relations and entropy current corrections 
in presence of minimal torsion, and compute constraints on hydrodynamic transport. In 
the following we only consider parity-odd sector; in parity-even sector we did not find any 
discrepancy between entropy current and equilibrium partition function to start with, and 
moreover calculation with torsion turns out to be trivially equivalent to what was done without 
torsion. 

Odd Dimensions (d=2n-l): In odd dimensions introduction of X^ does not lead to any 
new data; hence constitutive relations eqn. (3.14) do not modify, neither does the entropy 
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current correction eqn. (3.28). However divergence of entropy current does modify, 


Y7 _ iM 

— V M I c — IfcwTi 


M^s — ''( 0 ) 
n—1 


— 9 mWs( 0) + \^M — — C' 2 ,(l) {^MNU^ “ 

—dM^s{r) + ^s{r)^dM'& + f^C'2^(,.) — C'2,(r+1)'^ {^MNU^ — 

r=l L ^ 


+ I ^s{r) + 


E + P — 'dwR 


C'2,(r) — '&^C2,{r+l) ) ( 


(A.3) 


Plugging in the canonical part of entropy current eqn. (A.2) and demanding > 0, we 

will find a consistency condition on entropy current, 

^Ws(r) = ^ = 0 ; ~ ^ ^ (A.4) 

where C*i,(r-) is a constant. This is the missing constraint, which we did not find through 
entropy current in absence of torsion. Using this we can check that we get all the constraints 
eqn. (3.19) which we got from equilibrium partition function. One can now take a torsionless 
limit, which is trivial as there is no Xm dependence in the constitutive relations or entropy 
current. 


Even Dimensions (d=2n): The even dimensional case is more interesting. Constitutive 
relations eqn. (3.15) modify as. 


n—1 


frMN _ 


E OcmH-) - E 


r=0 


r=0 


n—1 


cM _ (n-l\iMN 

^{n) - I r )\r) 

r=0 
n—1 


Kvo{r)dN'^ + Ku{r,s)dNT^ + i^e{r)dN'& + Kx(r)^N + ^e{r) “ '^Onu) 


^ ^ (”r \vj{r)dN'^ + \u{r)dNV + kq{r)dN'& + \x{r)XN + ^q{r) {J^NRU^ “ 'OdNv) 

(A.5) 


r=0 


On the other hand, most generic entropy current correction (omitting terms that will give 
pure derivative terms in the divergence) will be given as. 


rf = E 


(r) 


uKT AD 


r=0 


n—1 


+ tfnE{V)'")'' 

r=0 




(A.6) 
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Its divergence can be computed to be, 




+ n ^ (%9MS2^{r) - dMS2,{r+l) \ {^NRU^ - 

r=0 L V / 

+ (^I(r)M'^ H-- dMS2^{r) — 'dl^(9M‘52,(r+l)^ 


E + P — -dwR, 


f d ~ d ~ \ 

^ v d'O D'uj 


“ 1 “^ ( p\ '^5'i9(r) o Q^si'{r) 1 d]\4ud]\[l^ H- I? ^ '^s^'(r) o '^5^A7(r) 1 ; (-^•^) 


where, 


^{r) ^'^5X(r) ~l~ ^szu(r)^'^ ^si'(r)^^ 


sv'{r)^^ “T ^S'i9(r)^ 


The divergence of canonical piece on the other hand is given as: 


79X7 

^}LM^s{can) 


-eEOCMt 


y~l ^fw{r)dM'^ + Kv{r)dMR + {^e{r) + 


+ i^tir) — Kx{r) ) {J^MRU — 'dduv) ( Xn - a^N'& 


+ y\-w(r)dM'^ + \v(r)dMR + ^^q(r) + ;^\x(r)J dM'&j {XnRU — I^Onr) ■ (A.9) 

Comparing the two we will get three consistency conditions on the entropy current, 

d ~ d ~ d ~ d ~ d ~ d ~ 

-^Xuj{r) = -Q;^^s'd{r)7 -^^s^{r) = -^Ku{r)7 '^^suir) = ^Vro(r)) (A.10) 

which will have most generic solution, 

Q 

^sro(r) fl(r) ; 

d d 

Ku{r) = + T;-f2{r){^,l^) 


■'suyri Q^J j - T- j ' j ^yr) y- / 

d d d 

+ ^f2{r){^,^) + ^/3(r)(^), 
for some arbitrary functions ro, i/), /2(,,)(ti7, i/), fs(^r){^)- We define. 


(A.ll) 


‘5l,(r) — + fl{r) + /2(r) + /3(r)- 


(A.12) 
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Table 3: Leading Derivative Order Data for Galilean Fluid (Non-Covariant) 


Newton-Cartan Data 

Non-Covariant Data 

Parity Even 


'Wiu'- - ^(e*^ - aiU^)gijdtg"^ - aidtu^ 

p^'^du'd, p^^'dyV 

- a* did, - a^dtw, VV - a^dtu 

pf^’^ {F^puP - 

e* + - dVV + §a^dtv 


-2a^^dtu^^ + ^g^^aidtu^ 
-(e'*’ - Oju*) {dtg'-^ + gkAg’"^) 

Parity Odd — Odd Dimensions (d = 2n — 1) 



Parity Odd — Even Dimensions (d = 2n) 

1 

< 

3 

< 

< 

* 

III 

o 



^(r) ~ ^(r) 

^\r) - (PidF) 

ffT) {F^pv!' - dd,u) 

^(r) (®i + ~ + ^ajdtiy) 


9kl 

where. 



One can check that = d5i (r). Demanding entropy current divergence to be non-negative 
one can check that (in torsionless limit) we get all the equilibrium partition function con¬ 
straints, as well as the additional entropy current constraints which we found before. 

Hence we have established the agreement of equilibrium partition function and entropy current 
in arbitrary number of dimensions. These results upon reduction, agree with the entropy 
current calculation for 2 spatial dimensional fluid in [2], except for a few discrepancies. We 
provide a detailed comparison with their results in appendix (C). 

B Non-Covariant Resnlts 

We mentioned in the main text that vector field T defines a reference frame. We can go 
to local rest of one such frame by choosing a basis = {x~,t,x^} such that V = d- and 
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T = dt- This essentially amounts to setting n* = 0 in the Newton-Cartan construction. For 
example in the equilibrium configuration discussed in § 2.2, we have studied the system in 
local rest of frame defined by isometry . Using the same field decomposition as given in 
§ 2.2, partition function eqn. (4.16) reduces to: 


6W = rI {dx^} - e'^^ (j* - - ^e'^A) 6ai + 


+ {p6Bt + + {qSAt + e-'^4<5A) 


(B.l) 


It is exactly same as eqn. (4.19), except that the observables are now defined with respect to 
frame T and are not independent of t. Note that choosing a frame fixes the Milne invariance 
in the partition function, and hence Milne Ward identity p* = goes on-shell. Under Milne 
boost background fields transform as, 




Bi ^ Bi 'ipi. 

On the other hand various densities and currents transform as. 


(B.2) 


P-^P, fp^fp-p^\ ^ 


1 - 


1 


etot etot - fpi^i + U “ UoA*) - A + (jp - A*) • 


ypv. . 2 ^ 

q^q, fq^jq-qi>'- 
The conservation equations in non-covariant basis becomes, 

-^dt (Vpgnc) + Vi (e“'^j*) = 0 , 

v9 


(B.3) 


V9 


dt iVgpnc) + Vi (e Ap) = 0 


■^dt {^/getot,nc) + Vi (e + e - e '*’a(T)ijp 

-^dt i^/gpnc,i) + Vj = -^e~'^aiB’^dtgjk + e"'^ {qet -b -b e""^ (-pa(T)i + ^(T)ijjp) 


where non-covariant densities gets a contribution from temporal curvature, 


(B.4) 


Pnc — P 6 jpCtii Uo4,nc — 6 j^CLii q-ac — q 6 jq^ii Pnc — Jp ® t ^CLj. (B.5) 

These are just the usual Galilean conservation equations, generalized to curved space-time. 
Constitutive relations of a Galilean fluid written in mass frame can be found as, 


ji = p* = Ru\ td = RvPu^ + Pgd + Ttd ^ ji = QpP + 
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(B. 6 ) 


Ru^j + 4 + TT^^Uj. 

s = S-T,_, jl = Su^ + ^4 - u4 + Ti - Ts-u\ 

They follow conservation laws eqn. (B.4). Finally we can explicitly obtain the constitutive 
relations upto leading order in derivative expansion for by reducing results of § 5 in mentioned 
basis. Reduction of various data down to non-covariant basis is given in table (3). In the 
following we present results for a special case when time is flat (a* = <I> = 0 ), space is time 
independent {dtgij = 0 ) and reference frame is inertial a* = ojij = 0 for simplicity. 


1 _9 

Etot = E + , 


ji = [E+p+ 


Odd Spatial Dimensions: We hrst present constitutive relations for a fluid living in odd 
spatial dimensions {d = 2n — 1), 

vr*-^ = — g'^^C'WkU^, 

n—1 

4 = + dKg (E + ^ 

r =0 

n—1 

4 = + ag {E + ^ Cb (B.7) 

In the special case of 3 spatial dimensions we will get the well known results [3], 

TT^d = -rja^d _ (g^d^kU^, 

4 = + a, (e* + (h X By - t9VV) + 

4 = Kg^^d + ag {4 + (u X By - WV) + + u;g(^yB\ (B. 8 ) 

where R* = ^ By = e^'d^VjBk = fd'^dy.^ 


Even Spatial Dimensions: Similarly in even spatial dimensions {d = 2n) we can get, 

n—1 / n \ 

-9« CV,s‘ + Oc,r,;,„ . 

r=0 V r=0 / 


r=0 
n—1 


r=0 




r=0 
n-1 


E -d P — dwR 


R 


djS2Xr) - ^k'djS2,y.+i) I , 


n-1 


?(J — ( 5 ^kiq + ^ ^ )^^r)^q{r) j djd + ( g ^Gg + ^ ^ )^(l)^q{r) j (^j + PjkU ddjV 


r=0 
n—1 


r=0 




(B.9) 
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These results for sure are messy, but take a cleaner form in 2 spatial dimensions, 

TT*^ = -ria^J - Qki - [c'Wku’" + + CbB^ , 

= [g'^l^e + dj'd + Kq - Kq) {cj + ejku’" B - 

+ ?9e*^ ^^^ djS2X0) - '&’^djS2Xi)^ , 

gq = {g'^Kq + e*-^kg) dj'd + {g'^dq + e*-^dq) (ej + ejku'^B - -ddju^ 

+ 52^(0) - a,52,(1)) , (B.IO) 

where B = Pij, ui = e^^diVj, e^^VjB = p'^^Vj. 

Under the assumption of flat time, time independent space and reference frame being inertial, 
these constitutive relations follow very simplified and familiar conservation laws, 

dtQ + '^ifq = 0, dtR + ^i{RvP = 0, 

dt {e + ^Rv\^ + ^di = j\ei, dtiRd) + = QE - jqjP^P (B.ll) 

Last term can be seen as (jq x By or e'^^jqjB depending on the number of dimensions. 

C Comparison with Geracie et al. [2] 

Our null fluid construction is computationally similar to [2], but has a different essence to 
it. Authors in [2] considered an extended (d + 2)-dim representation of the Galilean group, 
and realized it with the help of a (d + 2)-dim flat space. They then defined an extended 
Vielbein which connects this (d + 2)-dim space to (d + l)-dim Newton-Cartan manifold, and 
used this formalism to write Ward identities and constitutive relations of the Galilean fluid in 
a covariant manner. In this work however, we do hydrodynamics on the (d + 2)-dim curved 
manifold (null background) to start with, and later perform light cone reduction to get a 
Galilean fluid. As we mentioned, fluid on null background (null fluid) is essentially equivalent 
to the Galilean fluid, so we can expect computational similarities with the construction of 
[ 2 ]. 

[2] also studied torsional Galilean fluid in 2 spatial dimensions. We should be able to re¬ 
produce their results restricted to torsion-less case. In two spatial dimensions (d = 2), the 
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hydrodynamic frame invariants in parity-odd sector eqn. (5.14) become, 


(c 0 + Leo + /bB) , 

= ^arp^'^ {T^pvP - dd^v) + KTP^'^d^^ - {B^puP - ^d^v) + 

— fhtP'^B^puP + ^ 

li 

[BypuP - 'ddyu) -h aEpP'^dy'd + aE^P’^ {BypvP - ddyu) -t- eP'^ardy^ 

+ eP^dyrh + %eP‘'dyn, (C.l) 

H 

where we have made some redefinitions to make results look similar to [2], 


1 

II 

^ — *52,(0) 5 

d'T 

1 

= Kg- -m, 

rhe = 19^51(0) + 

2 

KT = Ke- -me, 

V 


aE = 

d(j. 

Cr* 

b 

II 

b 

ax' 

(C.2) 


~ - [^2^^ d dP d ^ dP dl 

U-Cu:--^^ 

7 7 \ n2 9P 9 dP 9 9P 5 1m 

[ dEd^ 9R9w 9Q9u\ ^ 


(C.3) 


Relations eqn. (C.l) are same as eqn. (6.64) of [2], except few subtle points. We do not have 
the coefficient kq while they don’t have h and f^j. We would like to mention few points which 
might explain this discrepancy. In the following we use notation as used in their paper. 


Just before eqn. (6.52) of [2], authors dropped the CuiLouP term in the entropy current as it 
gives rise to a ‘genuine second order data’ ui. However one can show that just like magnetic 
field, there exists an independent combination. 


Co; {OJUP + Cty) , 


(C.4) 


which has composite divergence. Here = uP\/ pvy and VP = eP'^Vy is the duality operation. 
In other words u = uP9pUJ is not an independent genuine data, and can be decomposed using 
first order Ward identities. 


cu = -a ;0 - eP^VpUy 

= (^lB-u?je- eP'^VpCy + ^^B-VP Ep - VP Gp - ^VPpVyp. (C.5) 

Therefore [2] missed the n and its dependent coefficients. For other discrepancy we note 
that they have a term —^{ac — kE)EPGp in entropy production eqn. (6.55), which implies 
dc = He- This will give rise to another consistency condition, which can be read out directly 
from eqn. (6.62), 

CQ + T {dQCT - 9tcq) = 0 => CQ = T[dQfi{T,EQ) + dQf2{i'Q)], ct = T5r/i(F, z^g), 

(C.6) 
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for some arbitrary functions /i(T, j/q), / 2 (t'Q)- One can in turn absorb these functions by 
redefining, 

(C.7) 

As an effect, in eqn. (6.62) of [2] all the c’s drop out. Consequently by eqn. (6.70) all c 
dependent coefficients kq, vanish (or can be absorbed in definition of m^). Barring 
these modifications, we find our results to be in exact correspondence with [2] for torsionless 
fluids. 


D Conventions of Differential Forms 


In this appendix we will recollect some results about differential forms, and will set notations 
and conventions used throughout this work. An m-rank differential form on A4(rf+2)) 
can be written in a coordinate basis as. 


1 


At 


A dx'^^ A ... A 


ml 


(D.l) 


where fi is a completely antisymmetric tensor. On Mi^d+ 2 )^ volume element is given by a full 
rank form, 


,{d+2) _ 


7eMiM2...Ma+2dx^^ A dx^^^ A ... A dx^‘^+^, 


(^ _|_ 2)!.. .. .... ’ 

where e is the totally antisymmetric Levi-Civita symbol with value eo,i, 2 ,...,(i+i = \/|C| and 
G = det Gmn- Using it, Hodge dual is defined to be a map from m-rank differential forms to 
(d + 2 — m)-rank differential forms, 

* = {d + 2- m)\ ) dx^^ A... A dx^'*+2—. (D.3) 

One can check that ** = sgn(G)(—. For us obviously sgn(G) = —1 due to 

Minkowski signature of the metric, but we tag along this factor for clarity. The exterior 
product of a differential form is defined to be, 

1 




{m + 1)! 

Integration of a full rank form is defined as, 


[(m -A 1)9 [Mi/^M2...m^+i]] dx^i A ... A dx^^+b 


(D.4) 


[ /x('='+2) = sgn(G) f {dx^} ^/]G\ * 

J^{d+2) J 

= sgn{G) j {dx^]^\ 




(d-|-2)!" r‘ivii...ivid+2- (U-5) 

Here the raised Levi-Civita symbol has value g 0 .i, 2 ,...,(i+i _ ggn(G)/Y^|G|. Integration of an 
exact full rank form is given by integration on the boundary. 


^A^(d+2) 9M{d+2) 


(D.6) 
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where given a unit vector N normal to boundary, volume element on the boundary is defined 
as = *N. 


Newton-Cartan Differential Forms 


A Newton-Cartan differential form is a differential form on ■M.{^d+ 2 ) which does not have a leg 
along V, i.e. Such a form can be expanded as, 

A dx^2 a ... A dx^-. (D.7) 

Volume element of NC manifold is defined as, 

^(d+i) (d + 1)! dx^a+i_ (D_8) 

Since there is no degenerate metric on NC manifold we cannot define a Hodge dual. Hodge 
dual can however be defined if we chose a frame T. We can hence define spatial differential 
forms with the requirement that they should not have any leg along V and H(r)- For these 
forms, indices can be raised and lowered using p/^‘' and We can define a spatial volume 
element, 

e('^) = A F] = i dx^^ A ... A dx^^ (D.9) 

and corresponding to it a Hodge duality operation. 




(m) 


— ic 


V Au A fj, 


(m) 


{d — m)\ \m\ 

One can check that ** = —sgn(G)(— 






dx"^ A ... A dx'^'^-"^. 

(D.IO) 


Non-Covariant Differential Forms 


Going to the local rest of frame T used to define the Newton-Cartan spatial forms, we can 
check that the spatial forms behave covariantly on the spatial slice, i.e. can be expressed 
as, 

^ ^ ^ (D ll) 

m\ 

One can check that the volume element defined before is indeed a full rank form on the 
spatial slice and can be written in this setting as, 

^ {V^V^eMm,...i,) dx'i A ... A dx'C (D.12) 

The Hodge dual * associated with it serves as Hodge dual operation on the spatial slice. 




(m) 


{d — m)\ \m\ 




ll.-.l-mj! •••Jd — 'n 


dx-^^ A ... A dx-^'*"’". 


(D.13) 
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Finally a full rank spatial form can be integrated on a spatial slice, 


J M(d) 


sgn(G) 


'M 


"V ^V ^ = sgn(5r) / {dx^} ^J\g\ * 


(<i+2) 


(D.14) 


Here g = det < 7 *^ = e^^y = —e^^G. Since gij is a spatial metric sgn(g') = +1. Other 
conventions and notations are same as relativistic case. 
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